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SUPERSONIC  FLUTTER  OF  SIMPLY 
SUPPORTED  ISOTROPIC  SANDWICH  PANELS* 

By  Larry  L.  Erickson  and  Melvin  S.  Anderson 
Langley  Research  Center 


SUMMARY 

fA  theoretical  solution  using  two-dimensional  static  aerodynamics  is  presented  for 
the  supersonic  flutter  characteristics  of  flat  rectangular  isotropic  sandwich  panels  with 
simply  supported  edges.  Tables  and  charts  giving  the  values  of  the  dynamic-pressure 
parameter  required  for  flutter  are  presented  for  various  values  of  panel  length-width 
ratio,  shear  flexibility,  and  midplane  stress.  It  is  found  that  a  decrease  in  transverse 
shear  stiffness  will  usually  lower  the  dynamic  pressure  required  to  induce  flutter.  How 
ever,  for  certain  cases  of  midplane  tension  the  opposite  effect  occurs.  Panel  mode 
shapes  are  also  presented  and  a  comparison  is  made  of  the  two-mode  Galerkin,  the  pre¬ 
flutter,  and  the  exact  flutter  solutions'^ 


INTRODUCTION 


Panel  flutter  is  an  important  design  consideration  for  vehicles  traveling  at  high 
Mach  numbers.  Consequently,  considerable  literature  has  been  published  dealing  with 
several  aspects  of  the  problem.  (See  ref.  1.)  One  aspect  which  has  not  been  adequately 
evaluated  is  the  effect  of  transverse  shear  flexibility  on  the  flutter  characteristics  of 
sandwich  panels.  Light-weight  structural  configurations  for  supersonic  and  hypersonic 
vehicles  may  incorporate  stiffened  panels  such  as  honeycomb  sandwich  panels  which  in 
most  cases  cannot  be  considered  rigid  in  shear.  Therefore,  unconservative  designs  may 
result  if  the  transverse  shear  stiffness  of  such  panels  is  not  taken  into  account. 

An  estimate  of  the  influence  of  shear  flexibility  on  panel  flutter  was  provided  in 
reference  2  for  flat  and  curved  isotropic  panels  and  in  reference  3  for  orthotropic  panels 
where  results  of  two-mode  Galerkin  solutions  were  presented.  However,  it  is  known  that 
these  approximate  solutions  become  increasingly  in  error  as  the  length-width  ratio 
increases.  As  in  references  2  and  3,  it  is  assumed  herein  that  the  aerodynamic  loading 
is  given  by  two-dimensional  static  aerodynamics  which  are  incorporated  with  the 

*The  basic  theoretical  development  presented  herein  was  given  in  a  thesis  by 
Melvin  S.  Anderson  in  partial  fulfillment  of  the  requirements  for  the  degree  of  Doctor  of 
Philosophy  in  Engineering  Mechanics,  Virginia  Polytechnic  Institute,  Blacksburg,  Virginia, 
June  1965. 


small-deflection  theory  for  flat  sandwich  panels  developed  in  reference  4.  However,  in 
the  present  investigation,  the  resulting  differential  equations  are  solved  exactly. 

The  results  of  reference  5  indicate  that  for  the  range  of  parameters  shown  therein, 
panel  flutter  analyses  based  on  this  simple  aerodynamic  theory  are  reasonably  accurate 
for  isotropic  panels  without  shear  flexibility  when  compared  with  the  results  obtained  by 
using  three-dimensional  unsteady  aerodynamics.  Thus,  for  most  cases,  the  results 
obtained  herein  are  not  expected  to  differ  greatly  from  results  that  would  be  obtained  by 
using  more  accurate  aerodynamic  theories. 

The  numerical  results  of  the  analysis  are  presented  in  the  form  of  tables  and  charts 
and  are  discussed.  Details  of  the  analysis  appear  in  the  appendixes. 


SYMBOLS 

A,B,C  coefficients  appearing  in  equations  (A6) 

A  =  - — - -  [kx  -  2n3  +  r(n3kx  +  (p) 

(l-rkx)Lx  ' 

a  length  of  panel 


B=rTTizB1  +  n2r)'n4] 


width  of  panel 


D 


D 


Q 


flexural  stiffness  of  isotropic  sandwich  panel, 

Eftthc2(i+Af 


Eftf3 

2(1  -  m2)  +  e(i  -  n2) 


,  tfV 

transverse  shear  stiffness  of  isotropic  sandwich  panel,  Gchc(l  +  —  j 

hcJ 


Dj^jDg  coefficients  defined  by  equations  (A15) 


Ej  Young's  modulus  for  faces  of  isotropic  sandwich  panel 

F(  )  function  defined  by  equation  (A25) 

Gc  shear  modulus  for  core  of  isotropic  sandwich  panel 


2 


depth  of  sandwich  core 


Lm,n 


integers 


midplane  stress  parameters 


7 r^D 


Z  lateral  aerodynamic  loading 

M  free-stream  Mach  number 

Mx,My  intensity  of  internal  bending  moments  acting  upon  a  cross  section  originally 
parallel  to  the  yz  and  xz  planes,  respectively 

Mxy  intensity  of  internal  twisting  moment  acting  in  a  cross  section  originally 

parallel  to  yz  plane  or  xz  plane 

in  exponent  in  equations  (A6) ,  denotes  roots  of  equation  (A9)  when  used  with 

subscript  j 

Nx,Ny  intensity  of  middle  plane  forces  parallel  to  x  and  y  axes,  respectively 
(positive  in  compression) 


Qx,Qy 


intensity  of  internal  shears  acting  in  z-direction  in  cross  sections  originally 
parallel  to  yz  and  xz  planes,  respectively 


1  +  r 


1  2 

free-stream  dynamic  pressure,  -PaV 


shear  flexibility  parameters, 


7j  D  7T^P 
b2DQ  a2DQ 


coefficients  defined  by  equations  (A30) 


„  A 
s  =  P  *  -s 

7)* 


3 


k 


thickness  of  sandwich  face  plates 


V  free -stream  velocity  of  airflow 

w  deflection  of  middle  surface  of  plate,  measured  in  z-direction 

x,y,z  orthogonal  coordinates  (see  fig.  1) 

0!,6,e  assumed  components  of  roots  of  equation  (A9)  (See  eqs.  (All)) 


0  =  \/m2  -  1 


-  mi 

=  — 

7777 


coefficient  defined  by  equations  (B6) 


y  = 


Xpr 


4ir2(l  -  rk^ 

?  =  y  -  « 


V 

X,X' 

X  =Xr? 


panel  length-width  ratio,  - 


dynamic  pressure  parameters, 


2qb2  2qa^ 
0D  ’  |3D 


3  1  +  nfr 
1  -  rkv 


Poisson's  ratio  for  sandwich  panel,  defined  in  terms  of  curvatures 


i-  =  y  +  a 

Pa 

Pm 


free-stream  mass  density  of  air 


mass  density  per  unit  area  of  panel 


frequency  parameters,  (7^-)  +  n2  and  +  n2  —— — 

W  n2D  \wr  /  tt2D 
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column  matrix 

Subscripts: 

cl 

°o  evaluated  at  —  =  00 

b 

cr  denotes  flutter  value  of  parameters 

p  denotes  preflutter  value 

A  comma  followed  by  a  subscript  denotes  differentiation  with  respect  to  the 
subscript. 

THEORY  AND  ASSUMPTIONS 

The  configuration  analyzed  is  shown  in  figure  1.  It  consists  of  a  flat  rectangular 
sandwich  panel  mounted  on  simple  supports.  The  panel  core  and  face  materials  are 
isotropic;  hence,  the  panel  itself  is  referred  to  as  being  isotropic.  The  panel  has  a 
length  a  and  a  width  b  and  is  subjected  to  uniform  midplane  force  intensities  Nx 
and  Ny  (positive  in  compression).  The  supersonic  flow  at  Mach  number  M  is  over 
the  top  surface  of  the  panel  and  is  parallel  to  the  x-axis. 

The  analysis  of  the  configuration  is  based  on  the  small  deflection  theory  for  flat 
sandwich  panels  developed  in  reference  4.  This  theory  incorporates  the  effect  of  shear 
deformations  by  expressing  the  total  panel  curvature  in  the  x-  or  y-direction  and  the 
twisting  distortion  as  the  sum  of  the  contributions  made  by  each  of  the  internal 
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shears  (Qx,  Qy)  and  moments  (Mx,  My,  MXy).  The  resulting  force-distortion  equations 
can  be  solved  for  the  three  separate  moments.  Substitution  of  these  expressions  for  the 
moments  into  the  equations  for  equilibrium  of  moments  about  the  x-  and  y-axes  and  the 
equation  for  equilibrium  of  vertical  forces  yields  three  independent  equations  relating  the 
lateral  displacement  w  and  the  two  average  shear  angles  Qx/Dq  and  Qy/Dq.  These 
equations  include  the  in-plane  force  intensities  and  the  lateral  loading.  In  the  analysis 
presented  herein,  the  lateral  loading  is  comprised  of  the  inertia  force  and  the  pressure 
due  to  supersonic  flow  (given  by  two-dimensional  static  aerodynamics).  In-plane  and 
rotary  inertia  loadings  are  not  considered. 

The  boundary  conditions  imposed  are  those  of  simple  supports,  but  inclusion  of 
shear  effects  requires  that  a  third  boundary  condition  be  satisfied  in  addition  to  the  two 
usual  conditions  of  zero  moment  and  middle-surface  deflection  along  the  panel  edges. 

This  third  boundary  condition  depends  on  the  assumption  made  for  the  panel  support.  If 
the  support  is  applied  only  to  the  middle  surface  of  the  panel,  the  boundary  condition  is 
that  MXy  must  vanish.  If  the  support  is  assumed  to  be  applied  over  the  entire  thickness 
of  the  panel  the  shear  angle  Qx/Dq  is  zero  along  an  edge  parallel  to  the  x-axis  because 
there  is  no  x-displacement  of  points  at  the  boundary.  Similarly,  the  shear  angle  Qy^Dq 
is  zero  along  an  edge  parallel  to  the  y-axis.  The  last  boundary  condition  (shear  angle 
of  0°)  is  usually  the  more  closely  approached  in  practice  and  is  the  one  used  herein.  The 
exact  solutions  to  the  differential  equations,  subject  to  the  stated  boundary  conditions, 
lead  to  a  transcendental  characteristic  equation  from  which  the  panel  frequencies  and 
mode  shapes  can  be  determined  as  a  function  of  midplane  loads  and  dynamic  pressure. 

For  the  type  of  aerodynamic -force  approximation  used  herein,  flutter  is  known  to 
occur  only  if  variations  of  airflow  or  panel  parameters  can  force  a  coalescence  of  two 
panel  frequencies  (ref.  6).  The  locus  of  such  points  forms  a  flutter  boundary  that  sepa¬ 
rates  a  region  of  stable  motion,  where  all  frequencies  are  real  and  distinct,  from  a 
region  of  dynamic  instability  where  at  least  one  pair  of  frequencies  are  complex  con¬ 
jugates.  Details  of  the  exact  solution  together  with  a  practical  procedure  for  obtaining 
numerical  flutter  results  are  presented  in  appendix  A. 

RESULTS  AND  DISCUSSION 

For  a  panel  of  given  length  or  width,  the  flutter  value  of  the  dynamic  pressure  q  is 
a  function  of  the  length-width  ratio  a/b,  the  bending  and  shear  stiffnesses  D  and  Dq, 
respectively,  and  the  in-plane  load  Nx.  For  the  presentation  of  results,  these  variables 
are  expressed  in  terms  of  nondimensional  parameters. 

Numerical  results  of  the  analysis  are  presented  in  tabular  form  in  tables  I  and  n. 
Table  I  is  for  a/b  <  1,  where  flutter  values  of  the  dynamic -pressure  parameter 
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,  _  2qa3 

Cr  ;3D 


and  the  frequency  parameter 


7T4D 


are  tabulated  for 


2 

various  values  of  the  shear -flexibility  parameter  r*  a  D--  and  the  stress  parameter 

N  2  ,  ** 
kx  =  — .  The  flutter  values  Xcr  were  obtained  by  plotting  frequency  loops  (X'  against 

7T2D 

<p',  see  fig.  2)  for  constant  values  of  r',  k^,  and  a/b.  At  the  peak  of  the  frequency 
loop  X'  =  \'cr  and  any  increase  in  X'  produces  flutter.  The  flutter  values  of  <p' 
associated  with  X^r  incorporate  Ny  and  can  be  used  to  determine  flutter  frequen¬ 
cies  u>cr.  As  in  the  case  for  panels  which  are  rigid  in  shear,  the  flutter  value  of  q  is 
not  affected  by  Ny.  Table  II  is  for  ^  =  1  and  to  keep  the  numerical  values  of  the  tabu¬ 
lated  results  within  reasonable  size,  the  parameters  have  been  redefined  in  terms  of  b. 


Thus,  in  table  n,  flutter  values  of  Xcr  = 


w2  + 


are  tabulated 


for  values  of  r  = 


tt2D 

b2DQ 


and 


2(lb3  m  _Pmh 

cr"“^"  ^cr  '  7t4D 

Nxbz 

ICj.  =  —A — .  Note  that  all  results  are  presented  in  terms  of 


n2 

N 

tt2D  Y 


IT2  D 


the  shorter  dimension  of  the  panel.  Tables  I  and  II  also  contain  values  of  a  parameter  a 
which  can  be  used  to  calculate  flutter  mode  shapes.  (See  appendix  A.)  In  table  n,  results 
are  not  presented  for  any  value  of  a/b  greater  than  that  required  to  produce  Xcr  =  0. 


The  results  in  tables  I  and  II  are  presented  graphically  in  figures  3  and  4.  To  keep 

these  figures  within  reasonable  size,  (x ' ) ^ /3  is  plotted  against  a/b  when  0  =  a/b  =  1 
1/3  '  cr 

and  XCj.  is  plotted  against  b/a  when  1  §  a/b  S  °°  for  different  values  of  the  shear- 
flexibility  parameters  r’,  r.  The  flutter  boundaries  appearing  in  each  individual  figure 
correspond  to  different  values  of  the  stress  parameters  k^,  kx.  Figure  4  is  for  stress- 
free  panels  with  various  values  of  r*  and  r.  It  shows  that  X’,r  is  essentially  indepen¬ 
dent  of  the  panel  width  for  a/b  g  1/4  and  that  Xcr  is  essentially  independent  of  length 
for  a/b  S  10.  Figure  4  also  illustrates  that  a  significant  reduction  in  Xcr  can  be 
caused  by  shear  flexibility  when  the  panel  carries  no  in-plane  loads  in  the  x-direction. 

It  should  be  noted  that  for  r  =  0  (Dq  =  °°),  the  panel  is  rigid  in  shear  and  for  this  case 
the  results  agree  with  those  of  reference  6. 


Effect  of  Shear  Stiffness  and  Stress  on  Flutter  Boundaries 

Figure  5  illustrates  the  effect  of  shear  stiffness  on  Xcr  for  a  square  panel. 

For  a  compressive  force  in  the  flow  direction  (kx  >  0),  the  theory  predicts  that  as  the 
panel  is  made  less  stiff  in  shear  (r  increasing),  Xcr  will  decrease.  The  same  result 
holds  for  the  stress-free  panel.  However,  when  the  panel  is  in  tension  (kx  <  0),  there  are 
cases  when  theory  predicts  that  a  panel  which  is  flexible  in  shear  (r  >  0)  will  flutter  at  a 
higher  value  of  Xcr  than  a  panel  of  the  same  geometry  which  is  rigid  in  shear  (r  =0). 
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For  example,  when  kjj  =  -4,  it  is  seen  that  after  a  slight  initial  decrease,  Acr  increases 
steadily  as  the  panel  becomes  more  flexible  in  shear.  At  r  =  2,  Acr  =  1230  which  is 
over  one-third  larger  than  the  r  =  0  value  of  Acr.  That  an  increase  in  panel  stiffness 
can  lower  Acr  is  an  unexpected  theoretical  result.  Whether  such  a  physical  phenomena 
does  actually  occur  or  whether  the  theory  employed  does  not  accurately  represent  the 
panel  behavior  under  certain  combinations  of  shear  stiffness  and  tensile  loading  is  not 
known.  Consequently,  such  results  should  be  regarded  cautiously  until  refinements  in  the 
theory  (for  example,  consideration  of  rotary  inertia)  or  experimental  evidence  either  con¬ 
firm  or  refute  this  anomalous  behavior. 

The  analysis  presented  in  appendix  A  is  based  on  the  criteria  that  flutter  occurs 
when  two  panel  frequencies  coalesce.  Any  circumstances  that  cause  two  in-vacuo  fre¬ 
quencies  to  coincide  will  then  produce,  by  this  definition,  flutter  even  though  the  dynamic 
pressure  approaches  zero.  Inclusion  of  either  aerodynamic  or  structural  damping  in 
theoretical  flutter  analyses  of  panels  which  are  rigid  in  shear  has  been  shown  to  remove 
these  zero-dynamic-pressure  flutter  points  (ref.  5);  the  resulting  theoretical  flutter 
boundaries,  however,  still  do  not  compare  well  with  experimental  boundaries.  Thus,  to 
gain  an  idea  of  conditions  for  which  the  theoretical  results  can  no  longer  be  considered 
reliable,  it  is  important  to  know  the  combinations  of  r,  kx,  and  a/b  which  produce 
zero-dynamic-pressure  flutter  points.  These  combinations  can  be  determined  from  the 
in-vacuo  vibration  characteristics  of  the  panel  as  described  in  appendix  C.  The  primary 

results  are  shown  in  figure  6.  When  kx  =  —  —  — Acr  approaches  zero  as  a/b 

(1  +  r)2 

approaches  infinity.  If  kx  is  greater  than  -  —  ■  ■,  then  Acr  goes  to  zero  at  values 

(1  +  r)2 

of  a/b  given  by  equation  (C5)  in  appendix  C.  If  kx  is  less  than  2--t-X—  then  Acr 

(1  +  r)2 

never  reaches  zero  but  has  a  finite  asymptotic  value  as  a/b  approaches  infinity. 

Comparison  of  Exact  Results  With  Approximate  Solutions 

Modal  solution.-  Development  of  the  exact  solution  provides  an  opportunity  for 
evaluating  the  approximate  two-mode  Galerkin  solution  of  reference  2.  A  comparison  of 
flutter  boundaries  is  shown  in  figure  7  for  stress-free  panels.  Note  that  the  two-mode 
solution  becomes  less  accurate  as  r  increases.  Also,  as  in  the  r  =  0  case,  it  becomes 
less  accurate  as  a/b  increases.  At  a/b  =  °°,  the  two-mode  solution  predicts  that 
ACr  =  0. 

An  explanation  for  the  growing  inaccuracy  of  the  two-mode  solution  with  increasing 
shear  flexibility  is  indicated  by  the  exact  flutter  mode  shapes  for  square  stress-free 
panels  shown  in  figure  8.  (These  shapes  were  calculated  (from  eq.  (A29))  by  using  the 
values  of  a  given  in  tables  I  and  n.)  As  the  panel  becomes  weaker  in  shear,  the  point  of 
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maximum  amplitude  moves  toward  the  trailing  edge  of  the  panel.  Thus,  it  would  be 
expected  that  the  representation  of  the  true  mode  shape  by  only  two  terms  in  a  sine  series, 
as  was  done  in  reference  2,  would  lead  to  increasingly  inaccurate  results  as  r  increases. 

Preflutter  solution. “In  reference  7,  a  solution  to  the  differential  equation  of  motion 
for  simply  supported  panels  that  are  rigid  in  shear  was  made  and  it  yielded  simple  alge¬ 
braic  expressions  for  X  and  <f>.  This  solution  corresponds  to  a  point  on  the  frequency 
loop  which  is  not  necessarily  at  the  peak;  thus,  a  value  of  X  is  given  which  is  less  than 
or  equal  to  Xcr,  hence  the  name  "preflutter. M  When  transverse  shear  effects  are  con¬ 
sidered,  it  is  found  that  a  simple  preflutter  solution  still  exists  (appendix  B),  even  though 
two  additional  differential  equations  and  an  additional  boundary  condition  must  be  satisfied. 

As  can  be  seen  from  figure  9,  which  is  for  stress-free  panels,  the  preflutter  solu¬ 
tion  gives  values  of  X  that  are  very  close  to  Xcr  when  a/b  is  sufficiently  large.  For 
small  values  of  a/b  the  preflutter  solution  is  often  in  poor  agreement  with  the  exact  so¬ 
lution.  Its  behavior  is  typical  of  that  shown  in  figure  10  which  is  for  r  =  2.0.  When 
a/b  =  1,  the  preflutter  solution  corresponds  to  a  point  low  on  the  second  leg  of  the  fre¬ 
quency  loop.  As  a/b  increases,  this  point  moves  along  the  second  leg  until  at  a/b  =  4 
it  is  nearly  at  the  top  of  the  frequency  loop  and  is  thus  virtually  identical  with  Xcr.  As 
a/b  continues  to  increase,  the  preflutter  point  moves  to  the  top  of  the  loop  (becoming 
exactly  equal  to  Xcr)  and  then  starts  down  the  first  leg.  Once  on  the  first  leg,  however, 
it  stays  relatively  close  to  the  top  of  the  loop,  and  thus  gives  a  close  approximation  to 
XCr* 

For  increasingly  large  values  of  a/b,  it  becomes  more  difficult  to  obtain  numerical 
results  from  the  exact  solution  for  Xcr.  Hence,  it  is  desirable  to  be  able  to  use  the  sim¬ 
pler  preflutter  solution  in  the  range  where  it  agrees  closely  with  the  exact  solution. 

Table  III  shows  a  comparison  of  the  preflutter  solution  with  Xcr  at  a/b  =  20  for  the 
values  of  r  and  kx  which  are  presented  in  table  II.  In  all  cases  the  solutions  differ 
by  less  than  2  percent.  Thus,  the  preflutter  results  used  in  table  H  and  in  figures  2  and  3 
for  a/b  >  20  are  justified  as  being  good  approximations  to  Xcr. 

CONCLUDING  REMARKS 

Charts  have  been  presented  to  facilitate  the  determination  of  theoretical  values  of 
the  dynamic  pressure  required  to  produce  flutter  of  flat  rectangular  isotropic  sandwich 
panels  with  simply  supported  edges.  The  flutter  value  of  the  dynamic  pressure  for  a 
panel  of  given  length  or  width  is  found  to  be  a  function  of  the  length-width  ratio,  the 
bending  and  shear  stiffnesses,  and  the  in-plane  load  acting  parallel  to  the  airflow.  It  is 
independent  of  the  in-plane  load  acting  perpendicular  to  the  airflow  as  long  as  buckling 
does  not  occur. 
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The  theory  predicts  that  a  reduction  in  transverse  shear  stiffness  usually  causes 
a  panel  to  be  more  susceptible  to  flutter.  An  unusual  result  is  obtained  for  panels  which 
carry  tension  loads  in  the  direction  parallel  to  the  airflow.  In  these  cases  it  is  theoret¬ 
ically  possible  for  a  reduction  in  shear  stiffness  to  make  a  panel  less  susceptible  to 
flutter.  However,  this  result  should  be  regarded  cautiously  since  at  present  there  is  no 
experimental  evidence  available  with  which  to  compare  such  theoretical  behavior. 

A  comparison  is  made  between  the  exact  solution  for  the  flutter  dynamic  pressure 
and  two  approximate  solutions.  A  two-mode  Galerkin  solution  is  shown  to  become 
increasingly  inaccurate  as  a  panel  becomes  more  flexible  in  shear.  A  relatively  simple 
preflutter  expression  is  found  to  give  results  that  are  in  good  agreement  with  the  exact 
solution  for  sufficiently  large  length-width  ratios. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  November  8,  1965. 
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EXACT  SOLUTION  FOR  THE  SUPERSONIC  FLUTTER 
BEHAVIOR  OF  SANDWICH  PANELS 

Differential  Equations 

The  small-deflection  equilibrium  equations  used  herein  to  represent  the  behavior  of 
a  sandwich  panel  are  derived  in  reference  4  by  use  of  the  following  force-distortion 
equations  which  relate  moments  Mx,  My,  MXy,  shears  Qx,  Qy,  and  displacement  w: 


x  8x\9x  Dq/  M  9y\ay  Dq 


%  -  °  ay  lay  Dq/+  M  dA  ax  Dq 


_  1  -  Md  8  /gw  Qy\  _9_Mv  Qx 
2  ax\8y  "  Dq/  8y^8x  Dq 


where  the  bending  and  shear  stiffnesses  as  given  in  reference  8  are 


Eftfhc  1  + 


Dq  =  Gch 


2(1  -  if  e(i  -  nf 

■44] 


respectively,  and  p  is  Poisson’s  ratio.  The  assumptions  involved  in  equations  (Al) 
require  that  a  straight  line  perpendicular  to  the  undeformed  middle  surface  of  the  panel 
remain  straight  and  of  constant  length  after  deformation  but  not  necessarily  perpendicular 
to  the  deformed  middle  surface.  This  inclination  in  the  x-  (or  y-)  direction  from  a  right 


DQ  d< 


angle  is  the  average  shear  angle  —  or  -f-  . 

dq  \  dq/ 

The  lateral  aerodynamic  pressure  given  by  static  linearized  two-dimensional  super¬ 
sonic  flow  theory  is 

(A3) 

where  q  =  -  paV^  is  the  free-steam  value  of  the  dynamic  pressure  and  /3  =  -  1. 

2 
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Substituting  the  above  expression  for  l ,  together  with  the  lateral  inertia  loading,  into  the 
three  equilibrium  equations  of  reference  4  yields 


-Nxw 


-  N„w  XTXT  +  Qx,X  +  Qy,y 


p  w 
Hm 


2q 


,tt  -  —  "  ,X 


w  Y=  0 


Qx  1 


-w,xyy  ”  w,xxx  -  —  + 


D  Dq 


Qx  + 
X.xx 


Qy  .  l 

,xxy  ,yyy  “  d  Dq 


Q, 


1  -  M 
2 

1  -  i) 


y,yy  \  2  /  y,xx 


Qy  + 

,yy 


Q, 


£ 


+ 


2 

l  +  p 


Qy 


2 


,xy 


xy 


=  0  >  (A4) 


=  0 


where  pm  is  the  mass  per  unit  area  of  the  panel  and  Nx  and  Ny  are  positive  in 
compression. 


Satisfaction  of  Boundary  Conditions  Along  the  Streamwise  Edges 

For  the  simply  supported  edges  parallel  to  the  x-axis  at  which  the  support  is 
applied  over  the  entire  thickness,  the  boundary  conditions  are  (see  ref.  4) 

w  =  0 


My=0 


Q, 


=  o 


'Q 


(A5) 


General  product  solutions  for  the  lateral  deflection  and  the  shears  which  satisfy 
these  boundary  conditions  are 


-x 

ma-,_  ..  y  Oicot 


w(x,y,t)  =  Ae  “"sin  nn  -  e 

b 


fig 


Qx(x,y,t)  =  Be  asin  nw  £  elwt 


(A6) 


Qy(x,y,t)  =  Ce  acos  m  ei^t  J 

where  n  is  an  integer  indicating  the  number  of  sinusoidal  half-waves  in  the  y-direction 
and  co  is  the  panel  frequency.  The  expressions  for  w,  Qx,  and  Qy  also  satisfy  the 
differential  equations  (A4)  provided  that 


12 


APPENDIX  A 


where 


Nontrivial  solutions  are  obtained  by  equating  the  determinant  of  the  matrix  in  equa¬ 
tion  (A7)  to  zero.  Expanding  this  determinant  leads  to 
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where 


y 


\7]T 


4tt2(1  -  rkx) 


A  =  ”2  r 


1  -  rkx 


jk^.  -  2n2  +  +  <pj 


-  .  q  1  +  n2r 

X  =  \r)6-. - r- 

1  -  rkY 


B  =  — 2 - L(l  +  n2r)  _  n4 

1  -  rkvL 


(A10) 


If  the  roots  of  equation  (A9)  are  written  in  a  form  similar  to  that  used  in  refer¬ 
ence  6, 


m^  =  y  +  a  +  id 


m2  =  y  +  a  -  i6 


m3  =  y  -  01+  e 


m4  =  y  -  o?  -  e 


mg  =  Tirjn 


1  + 


(1  -  m)  n‘ 


1/2 


m6  =  "m5 


J 


(All) 


they  can  be  related  to  the  coefficients  in  equation  (A9)  as  follows 


6y2  -  2c?  +  6^  -  =  t?A 


(A12a) 
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2a  (b^  +  e%)  -  2y  6^  -  +  2 (y2  -  aty 


=  X 


(A12b) 


y4  +  y2^j2  _  g2  _  2a^j  _  2ay(g2  +  5^j  +  (a?1  -  e^j{a^+  6^)  =  -ir^B  (A12c) 


For  computational  purposes  it  is  convenient  to  solve  equations  (A12a)  and  (A12b) 
for  6  and  e . 

•\ 

.27 


§2  =1.+  a2+  —  -  £  2y2  +  3ya  -  tt2£ 

4a  O  Gt\  y 


2=2._02.»?A.t(2  2_3  2A 

4a  2  <A  2, 


(A13) 


It  should  be  noted  that  62  and  could  be  negative  and  this  possibility  was  allowed 
for  in  the  calculations  by  taking  into  account  any  imaginary  quantities  that  appear.  How¬ 
ever,  for  the  ranges  of  parameters  covered  in  the  calculations,  6  was  real  for  any 
point  on  the  frequency  loop  and  e  was  real  except  at  sufficiently  high  axial  compression 
where  e  became  imaginary.  Using  these  relations  to  eliminate  6  and  e  from  equa¬ 
tion  (A12c)  yields  the  following  cubic  in  a  2 


a®  -  Dja^+  Dga^  -  Dg  =  0 


(A14) 


where 


»3 


(A15) 


Equations  (A13)  and  (A14)  provide  a  procedure  for  obtaining  numerical  values  of  X  once 
the  boundary  conditions  at  x  =  0  and  x  =  a  are  satisfied. 
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Satisfaction  of  Boundary  Conditions  Along  the  Leading  and  Trailing  Edges 

For  the  simply  supported  edges  parallel  to  the  y-axis  at  which  the  support  is 
applied  over  the  entire  thickness,  the  boundary  conditions  are  the  same  as  those  given  by 
equations  (A5)  if  y  and  x  are  interchanged.  Before  the  boundary  conditions  can  be 
applied  to  the  solutions  given  by  equations  (A6),  it  is  necessary  to  express  w,  Qx,  and 
Qy  in  terms  of  one  set  of  coefficients  (that  is,  either  A,  B,  or  C).  The  coefficient  A 
can  be  eliminated  from  the  second  and  third  rows  of  matrix  equation  (A7).  The  result  is 


n^r/l  -  p\/m j 


(A16) 


When  j  equals  5  or  6,  the  first  bracketed  term  is  equal  to  zero  (eqs.  (A9)  and  (All))  and 
no  information  is  obtained  about  the  relationship  between  Bj  and  Cj.  However,  when 
j  =  1,  2,  3,  4,  the  first  bracketed  term  is  not,  in  general,  equal  to  zero.  Therefore, 


-  177711 R- 

Cj  BJ 

mj 


0  =  1,2,  3,  4)  (A17) 


Replacing  Cj  with  Bj  in  the  third  row  of  equation  (A7)  gives 


2?j2  +  r  (71^2 172  -  mj^  ab2 


j7r2  ^2^772  -  mj2^  D 


Bj  (j  =  1,  2,  3,  4)  (A18) 


If  Cj  is  eliminated  from  the  first  and  second  of  equations  (A7),  it  is  found  that 


44 + w + -  -4  4 1  -  - + 


(A19) 


Again,  the  bracketed  term  before  Bj  is  zero  for  j  equal  to  5  or  6.  The  coefficient 
of  Aj  is  not,  in  general,  equal  to  zero  and  therefore 


A5  =  A6  =  0 


(A20) 


Thus,  for  j  equal  to  5  and  6,  the  first  of  equation  (A7)  gives 


mj 

Cj  -  m 7?®] 


(A21) 


16 


APPENDIX  A 


Use  of  these  relations  enables  equations  (A6)  to  be  written  in  terms  of  Bj  alone  as 
follows: 


Substituting  the  above  expressions  for  w,  Qx,  and  Qy  into  the  boundary  con¬ 
dition  equations  yields 


where 


(A23) 
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The  nontrivial  solution  for  the  deflected  shape  of  the  panel  is  obtained  by  equating  the 
determinant  of  this  matrix  to  zero.  By  manipulating  rows  and  columns,  additional  zeroes 
can  be  introduced  into  columns  5  and  6  of  this  determinant  so  that  it  becomes 


1 

1 

1 

1 

e^l 

m2 

e 

em3 

rr 14 
e 

mi2 

iS22 

—  2 
m3 

—  2 
m4 

=  0 

(A2/- 

m^1 

H22e”2 

.  -  2  53 
m3  e 

m4^e  ^ 

Note  that  equation  (A24)  is  independent  of  the  roots  in5  and  m0.  It  is  mathematically 
the  same  as  for  a  panel  which  is  rigid  in  shear  except  that  the  roots  mj  are  functions 
of  the  shear  stiffness.  (See  eqs.  (All),  (A10),  and  (A8).)  Expanding  equation  (A24)  and 
replacing  the  roots  mj  with  the  expressions  given  by  equations  (All)  yields. 


F(a,6,e,y)  = 


(q2  +  g 2)  +  4^  ^2  _  g2^  +  +  52  -  e 2) 

85e(o?  -  y2)(cosh  e  cos  5  -  cosh  2a)  =  0 


sin  5  sinh  e 


(A25) 


Equation  (A25)  reduces  to  equation  (9)  of  reference  6  when  y  equals  zero  (infinite  trans 
verse  shear  stiffness). 


Relations  (A10),  (A13),  (A14),  and  (A25)  are  8  independent  equations  involving  13 
quantities  (r,  kx,  7 7,  0,  n,  A,  A,  A,  B,  y,  a,  6,  e).  For  simple  supports  these  equations 
constitute  the  solution  to  equations  (A4)  once  5  of  the  13  quantities  are  specified. 

The  differences  between  the  results  of  the  analysis  presented  herein  and  the  results 
of  reference  6  are  due  to  the  parameter  r.  In  figure  11,  some  of  these  differences  are 
illustrated  by  showing  the  relation  between  the  frequency  loops  and  the  functions  A 
and  B  which  are  defined  by  equations  (A10).  Regardless  of  r,  A  reaches  the  value 
Acr  when  two  of  the  frequencies  of  the  panel  coalesce.  An  increase  in  A  above  Acr 
causes  cn2  to  become  complex.  Thus,  one  of  the  two  square  roots  of  co2  must  possess 
a  negative  imaginary  part  which  by  equation  (A6)  is  associated  with  divergent  motion  or 
flutter. 

When  r  is  zero  (infinite  shear  stiffness),  the  expressions  for  A  and  B  reduce 
to  the  definitions  in  reference  6.  In  this  case  A  is  not  a  function  of  the  frequency 
parameter  (p  and  depends  only  on  n,  77,  and  k^..  The  frequency  loops  (variation  of  A 
with  <p)  then  lie  in  planes  of  constant  A  and  Acr  is  a  function  of  A  only.  Hence, 
differing  panel  configurations  with  the  same  value  of  A  but  with  different  combinations 
of  77  and  kjj  will  all  theoretically  flutter  at  an  identical  value  of  Acr.  For  r  greater 
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than  zero  (finite  shear  stiffness),  A  is  no  longer  independent  of  (p.  As  a  result,  the 
frequency  loops  are  rotated  out  of  the  constant  A  planes  by  the  angle 

^tan-ilA?  (A26) 

1  +  r 

In  this  case  n,  r /,  and  kx  do  not  form  a  single  parameter  determining  Acr. 

Two  additional  effects  of  nonzero  r  can  be  seen  from  figure  11.  One  of  these  is 
that  the  frequency  loop  tends  to  bend  over  in  the  direction  of  increasing  frequency.  As 
r  increases,  this  effect  becomes  more  pronounced.  The  second  feature  is  that  for 
A  =  0  (no  airflow),  the  natural  frequencies  of  the  panel  are  reduced  in  magnitude  as  r 
is  increased. 

Procedure  for  Obtaining  Numerical  Values  From  Exact  Solution 

Since  the  flutter  value  of  A  is  sought  for  a  given  panel  configuration  and  midplane 
stress  condition,  it  is  desirable  to  specify  r,  1^,  and  r\.  If  n  and  the  frequency 
parameter  <p  are  then  prescribed,  the  remaining  eight  quantities  can  be  determined. 
This  process  enables  one  to  obtain  the  variation  of  A  with  (p  for  fixed  values  of  the 
length-width  ratio  and  stress  and  shear -flexibility  parameters. 

A  practical  procedure  for  obtaining  flutter  values  of  A  for  selected  values  of  r, 
kx,  T],  and  n  is  as  follows: 

(1)  Select  (f>,  a  good  choice  being  a  value  midway  between  the  first  and  second  in- 
vacuo  values  of  (p  (see  appendix  C). 

(2)  Calculate  A  and  B  from  equations  (A10). 

(3)  Make  an  initial  estimate  for  the  correct  value  of  A.  Calculate  corresponding 
values  of  A  and  y  from  equations  (A10). 

(4)  Solve  equation  (A14)  for  o?.  The  Newton-Raphson  technique  works  very  well 
since  only  one  real  root  of  equation  (A14)  was  found  to  exist  when  the  estimated  value 
of  A  was  close  to  the  correct  value. 

(5)  Calculate  5  and  e  from  equations  (A13). 

(6)  Use  equation  (A25)  to  calculate  F(a?,6,e  ,y),  a  nonzero  value  meaning  an 
incorrect  choice  of  A. 

(7)  Repeat  steps  3  to  6  until  a  value  of  A  is  found  that  differs  from  the  correct 
value  (that  is,  F(a',6,e,y)  =  (5)  by  only  an  acceptable  amount.  (An  allowable  error  in 

A  of  0.01  percent  was  used  in  calculating  the  results  presented  herein.) 
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(8)  Increase  cf>  and  repeat  the  process  until  the  point  =  0  is  obtained. 

When  the  numerical  results  presented  in  this  report  were  obtained,  the  increments  in  <p 
were  taken  small  enough  to  assure  that  the  values  of  Xcr  were  in  error  by  no  more 
than  1  percent.  In  most  cases  they  are  correct  to  0.1  percent.  In  all  cases,  calculations 
show  that  the  critical  flutter  solution  is  obtained  by  setting  n  =  1  and  determining  Acr 
from  the  coalescence  of  the  two  lowest  frequencies. 


Panel  Mode  Shape 


The  lateral  deflection  can  be  obtained  from  the  first  of  equations  (A6)  once  the 
coefficients  Aj  are  known.  These  coefficients  can  be  determined  from  equation  (A23). 
If,  in  equation  (A23),  n  is  set  equal  to  1  and  rows  5  and  6  are  multiplied  by  r,  and 
row  5  is  subtracted  from  row  3  and  row  6  is  subtracted  from  row  4,  the  coefficients  of 
B5  and  Bg  appearing  in  rows  3  and  4  can  be  made  equal  to  zero.  Then,  replacing  Bj 

i//j0jAj  for  j  =  1,  2,  3,  or  4,  where 


with 


b3 


and  recalling  that  Ag  =  Ag  =  0  yields 


1 

1 

1 

1 

rA\ 

Ai 

V 

e“l 

eS2 

e™3 

e™4 

A2 

I 

0 

*1 

^2 

^3 

^4 

< 

A3 

H 

o> 

mi 

^e  1 

So 

2 

mo 

^3e 

mi 

^4e  4 

0 

S.  > 

(A27) 


(A28) 


By  a  corrollary  to  Cramer's  rule  (ref.  9),  the  coefficients  Aj  can  now  be  determined 
within  an  arbitrary  constant.  Substituting  the  expressions  for  Aj  and  the  expressions 
for  mj  given  by  equations  (All)  into  equation  (A6)  and  ignoring  the  arbitrary  constant 
yields,  after  considerable  manipulation, 
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i* 

So 


rX 

J  o 


e  acos  -  e  acosh 

ci  cl 


^  =  j^Te^sinh  e  +  (^3  -  i£4)e^sin  6 
— 

+  Se^sinh  e  +  (^3  -  i^e^cos  5  +  i/^-6  -  i/^e^ +e 


1  J-X 

|e  asin  6  tt 


^Te^cosh  e  +  (1^3  +  i^Je^sin  6  +  e^(T  cos  6  -  S  sin  6)je^asinh  ejj  (A29) 


where 


£  =  y  +  a 
%  =  y  -  a 


^-1 


S  _  (1  +  r)7T*?74  +  rG2  +  52)^  -  tt2??2(1  +  2r)(^2  -  52) 

(1  +  r^iAri^  +r^^2  +  §2^  _  277^r7^r (1  +  r)(|^  -  6^) 

T  _ _ -47T2T?2g6 _ 

(1  +  r^Ti^rft  +  +  6^)  -  2u2?7^r(l  +  r)(^  -  6^) 


(A30) 


The  flutter  mode  shapes  for  w  presented  in  figure  8  were  obtained  from  equation  (A29) 
by  using  the  values  of  o^.r  in  tables  I  and  II. 


21 


APPENDIX  B 


PRE  FLUTTER  SOLUTION 

From  equation  (A25)  it  is  seen  that  the  transcendental  equation  F(q'  5,e,y)  =  0  is 
satisfied  identically  when  e  =  2a  and  6  =  2m7r  where  m  is  an  integer.  (This  same 
relationship  was  noticed  in  reference  7  for  panels  which  are  rigid  in  shear,  and  leads  to 
the  so-called  "preflutter"  solution.)  Substituting  these  expressions  for  e  and  8  into 
equations  (A12)  and  setting  m  =  1  yields  the  following  expression  for  X 


_X_ 

tj-3 


-  A+  6 


'4  -  A 
6 


8ff 

l7T 


(Bl) 


(Setting  m  =  1  restricts  the  preflutter  solution  to  a  point  on  the  first  frequency  loop.) 

If  the  panel  is  rigid  in  shear  (r  =  y  =  0),  the  following  preflutter  equations  for  X  and  B 
are  obtained 


(B2) 


where  A  =  T?2(kx  -  2)  is  now  independent  of  the  panel  frequency.  As  noted  in  reference  5, 
the  first  of  equations  (B2)  yields  values  of  X  that  are  very  close  to  Xcr  when  A  is 
negative.  (See  fig.  9.) 

For  r  *  0,  preflutter  values  of  X  are  not  so  easily  obtained  since  it  is  not 
possible  to  solve  for  X  as  an  explicit  function  of  r,  77,  and  kx.  However,  it  is 
possible  to  write  equation  (Bl)  in  terms  of  r,  77,  and  A  by  replacing  £  with 

_x/r/? 7_V  (See  eqs.  (A10).)  The  result  is 
4tt3\1  +  r/ 


( 


1  +  RA  -  R2 


i  +  ra-b^-^M 


12 


”3T”  (1  +  4R) 


4R)(lO  -  a)2  (4  -  a) 


(B3) 
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where 


R  = 


r  Ai2 

1  +  r 


Thus,  if  r,  r],  and  A  are  specified,  X  can  be  calculated.  (Note  that  A  is  a  function 
of  the  unknown  frequency  parameter  $.)  Since  X  is  known,  can  be  determined  from 
equations  (A10);  and  then  2?  and  B  can  be  obtained  from  equations  (A12).  The  stress 
parameter  kx  is  obtained  by  eliminating  <p  from  the  expressions  for  A  and  B 
given  by  equations  (A10): 


kx  = 


?]2|JLI  +  (A  -  rb) 

1  +  r  2 _ _ 

77^(1  +  r)  +  r(A  -  rb) 


(B4) 


Finally,  X  is  determined  from  the  third  of  equations  (A10).  It  should  be  noted  that  for 
fixed  r  and  ??,  k^.  cannot  be  specified  arbitrarily  since  it  is  determined  by  the  choice 
of  A  used  in  equation  (B3).  Thus,  trial-and-error  calculations  involving  different 
choices  of  A  are  necessary  to  obtain  X  for  the  kx  of  interest. 

The  preceding  preflutter  equations  can  be  used  for  infinitely  long  panels  by  letting 
T)  approach  °°.  This  condition  leads  to  the  following  expressions  for  X  and  kx: 


1  -  rkx 


(B5) 
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where 


4  .sst-isa!/ 4V«*pn 


J1 


r  = 


p  = 


8 


w 


(pit  _  1  f_A 


®  \?72, 


3/2 


jl/aV 

1 2L2J 


—2  \  (12  +  12s  -  s2)  ±  |^144  +  288s  +  120s2  +  M  s2  + 


2776tt6; 


1  +  r 


s  =P 


V‘ 


12P' 


J 


(B6) 


Thus,  by  specifying  r  and  — ,  X  and  k^.  can  be  calculated  for  rj  = 

??2  _ 

A 

An  interesting  solution  occurs  when  —  =  0.  K  the  negative  root  is  used  to  cal- 

/r\  V2 

culate  T,  r  =  — -  =0.  Hence, 

WJoo 


kx  = 


2  +  r 

(1  +  r): 


X  =  0 


(B7) 


That  is,  as  1^  is  increased  to  the  value 


2  +  r 


,  the  preflutter  value  of  X  for  an 


(1  +  r) 

infinitely  long  panel  drops  to  zero.  It  is  shown  in  appendix  C  that  Xcr  also  goes  to 

2  +  r 


zero  as  k^.  approaches 


(1  +  r)2 


when  rj  =  °°. 
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PANEL  BEHAVIOR  FOR  ZERO  DYNAMIC  PRESSURE 

The  natural  vibration  characteristics  of  the  panel  can  be  obtained  for  zero  dynamic 
pressure  by  setting  A  =  X  =  y  =  0.  Then  by  equation  (A14),  a  is  also  equal  to  zero. 
Equation  (A25)  reduces  to  (e 2  +  62)^sin  5  sinh  e  =  0  so  that 

6  =  m7r  (Cl) 

where  m  is  an  integer  designating  the  number  of  sinusoidal  half-waves  in  the  x-direction. 

o— 

From  equations  (A12)  and  (Cl),  it  is  readily  verified  that  =  - —  and 

m2 

B  =  m2(m2  -  A)  (C2) 


Equations  (C2)  and  (A10)  lead  to 


(CS) 


Use  of  the  expression  for  <p  as  given  by  equation  (A8)  enables  the  in-vacuo  natural  fre¬ 
quencies  to  be  expressed  as 


Equation  (C4)  can  be  used  to  calculate  the  in-vacuo  buckling  loads  by  equating  the  fre¬ 
quency  to  zero. 

In  the  analysis  presented  in  appendix  A,  it  is  assumed  that  flutter  occurs  when  two 
frequencies  of  the  panel  coalesce  or  become  equal.  For  certain  values  of  the  stress 
parameter  kx,  two  of  the  panel1  s  in-vacuo  natural  frequencies  become  equal.  This  con¬ 
dition  occurs  whenever  two  consecutive  mode  lines  cross  one  another  as  in  figure  11 
(that  is,  eq.  (C2)  is  satisfied  by  two  consecutive  integer  values  of  m).  The  numerical 
calculations  which  were  made  indicate  that  Acr  =0  at  these  points;  that  is,  the  frequency 
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loop  degenerates  to  a  point.  The  combinations  of  r,  kx,  and  77  which  produce  these 
theoretical  zero-dynamic-pressure  flutter  points  can  be  determined  by  using  equation  (C3) 
to  eliminate  <p  from  the  expression  for  B  given  by  the  last  of  equations  (A10).  Using 
the  in-vacuo  relation  between  B  and  A  (eq.  (C2))  then  leads  to 


2  +  r 

(1  +  r)^ 


-  rkx\  —1  -  rkx 

- A  +  rB - -  =  0 

1  +  r  /  1  +  r 


(C5) 


This  equation  gives  the  zero-dynamic -pressure  flutter  points  in  terms  of  r,  kx,  and 
?7.  The  values  of  A  and  B  to  be  used  are  those  for  which  equation  (C2)  is  satisfied 
by  two  consecutive  integer  values  of  m.  (The  first  zero-dynamic-pressure  flutter 
point  occurs  at  A  =  5.)  (See  fig.  6.) 
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TABLE  n.-  FLUTTER  SOLUTIONS  FOR  PANELS  WITH  LENGTH- WIDTH  RATIOS  GREATER  THAN  OR  EQUAL  TO  ONE  -  Continued 
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355 


EFLUTTER  AND 

- - cl 


.9 

1.859 


.7 

1.4 

1.528 


''■cr 


249.7 

136.4 

88.87 

48.68 

17.54 

.0429 


249.8 

135.6 

87.31 

46.37 

17.57 

.0371 


250.2 

133.9 

83.59 

40.64 

16.77 

1.183 

.0248 


AP 


249.3 

136.0 

88.58 

48.45 

17.39 

.0422 


249.4 

135.3 

87.02 

46.15 

17.43 

.0364 


249.8 

133.6 

83.31 

40.46 
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.0244 
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2 
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0 

.5 
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1.224 


4 
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^cr 
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79.84 
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.0053 

.0053 

255.1 
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129.5 

70.78 

70.51 

17.05 

16.97 

11.24 

11.18 

5.911 

5.880 

1.474 

1.462 

.0016 

.0016 

Figure  7.-  Comparison  of  flutter  boundaries  from  two-mode  Galerkin  solution  and  exact  solution,  k, 
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